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Abstract 


Study of the electronic properties namely conductivity, mobility and dominant 
transport mechanism across the grain boundaries, in case of boron doped 
polycrystalline-silicon films, was carried out in this work. The study shows that 
the average carrier concentration varies with the level of boron dopant 
concentration in the material. While at high dopant concentration most of the 
dopants contribute to average carriers, the average carrier concentration at low 
doping concentration is only a small fraction of the dopant concentration. This 
phenomenon is attributed to the presence of grain boundary traps in material 
.the study shows that at Boron dopant concentration of 1 0 14 cm' 3 , the free carrier 
concentration for 200A° and 1220A° grain size material is approximately 10 9 cm‘ 3 
The study also shows that at lower dopant concentration (<10 17 crrf 3 ) thermionic 
emission (TE) is the dominant mechanism behind carrier transport across the 
grain boundaries. However, both thermionic emission (TE) and thermionic field 
emission with scattering (TFES) contribute to carrier transport at high dopant 
concentration (10 17 -10 18 cm' 3 ). Further, at very high doping concentration, TFES 
is the dominant mechanism of transport across the grain boundaries. 
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Chapter 1 
Introduction 


Polysilicon has many important applications in integrated-circuit technology. 
Heavily doped polycrystalline -silicon films are commonly used as the gate 
electrodes and interconnections in silicon-gate MOS integrated circuits. Lightly 
doped films are frequently used for high value resistors in static memory circuits. 
Large gain polysilicon has the potential for the large volume production of low 
cost solar cells. In the past years, particular attention has been given to 
polycrystalline -silicon films and their application in microelectronics. 

Because of the wide applications of polycrystalline -silicon in integrated circuit 
technology, its electrical properties have been the focus of increasing attention 
by many researchers. These investigations have shown that electrical properties 
of polycrystalline -silicon are markedly different from those of single crystal 
silicon. At low dopant concentration the resistivity is several orders of magnitude 
higher than that of single crystal silicon and is not sensitive to change in the 
dopant concentration .At medium dopant concentration a small increase in 
dopant concentration results in a sharp drop in the resistivity. At high dopant 
concentrations the resistivity approaches that of single crystal silicon, although it 
always remains slightly higher. 

This high resistivity of lightly doped polycrystalline -silicon has been utilized in 
various applications such as high value load resistors in static RAM’s. 
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During the last two decades, a significant amount of work concerning the various 
properties of un doped or doped polysilicon films at different annealing 
temperature was done. Most of the work was oriented towards phosphorous - 
doped polysilicon films used as the gate material in metal oxide semiconductor 
(MOS) technology and to arsenic -doped polysilicon mainly used as the emitter 
contacts in bipolar technology. Polycrystalline -silicon thin-film transistors (poly-Si 
TFT) have been studied because of their potential application in high - 
performance active matrix liquid crystal displays (AMLCD). 

Much less attention has been given to boron-doped polysilicon especially at low 
dopant concentration, which is now commonly used in the state of the art silicon 
self -aligned bipolar technology as the extrinsic base contact. Boron doped 
polysilicon has drawn considerable attention for its potential application in 
polysilicon devices namely n-channel thin film transistors (TFTS). Therefore, a 
better knowledge of its conduction properties is necessary for proper design of 
the device. 

In the present work a model based on a combination of, carrier trapping and 
carrier Tunneling at grain boundary is used to explain the electrical conduction in 
B-doped polysilicon. To Study the dominant mechanism of transport at different 
doping levels analytical calculation is done first and then the results are 
compared with the experimental results known from various references. 
Calculated results accurately predict the conduction properties of polycrystalline - 
silicon films as a function of dopant concentration. Special interest is shown to 
calculate free carrier concentration at very low doping concentration (<10 16 cm‘ 3 ). 
The report is organized in remaining five chapters with the following contents: 
Chapter 2 summarizes some basic definitions and theories, which are used as 
tools in the following chapters. 

In Chapter 3 analytical calculations of the transport properties, resistivity and 
mobility have been done. 
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In Chapter 4 numerical calculations of grain boundary width and results obtained 
are shown and a comparison is made with the experimental results. 

In Chapter 5 analysis and discussion of the results are given. 

Chapter 6 gives the conclusions and further scope to improve the results. 


Chapter 2 

Models for Conduction in Polysilicon 


2.1 Polysilicon Grain & Grain Boundary 

In the single crystal silicon, dopant impurities diffuse in to the perfect lattice 
structure by interacting with point defects, such as silicon vacancies and 
interstitials. A polycrystalline material is composed of small crystal lattices joined 
together by grain boundaries. Inside each crystallite the atoms are arranged in a 
periodic manner so that it can be considered, as a small single crystal .The grain 
boundary is a complex structure, usually consisting of a few atomic layers of 
disordered atoms (Figure-la -1b). In case of polysilicon, two silicon crystals are 
connected through the grain boundaries alternatively. Like open crystal structure 
associated with the dangling bonds at dislocations allows rapid dopant migration 
along the dislocations, likewise in the polysilicon grain boundaries disorder ness 
provides high diffusivity paths along which dopant atoms can easily move. Even 
though the grain boundaries occupy only a small fraction of the sample volume, 
dopant migration along these path can markedly increase the over all dopant 
diffusion in the polysilicon. 

The dopant diffusion depends strongly on the crystal structure of the polysilicon 
as well as the detailed structure of the grain boundaries. Although the dopant 
atoms are not electrically active while they are in the grain boundary after 
diffusing along the grain boundary, they can move back into the grain and 
influence the electrical properties of the polysilicon. Grain boundary is a thin 
region of high diffusivity, normal to the sample surface between the two grains of 
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low diffusivity. Enhanced diffusion associated with the grain boundary can be 
described as a two -step process, the dopant atom first diffuses rapidly along the 
disordered structure of the grain boundary and then moves laterally into the 
grain, where it diffuses less rapidly. 



Grain Boundary J 



Grain Boundary 


Figure la: Grains and Grain Boundaries 



Figure 1b: Dangling bonds at the grain boundary 


2.2 Models for Conduction Mechanism in Polysilicon 

At low dopant concentrations, polycrystalline -silicon is characterized by a 
resistivity much higher than that of single crystal silicon, with low sensitivity to 
changes in dopant concentration [1] - [7]. At medium dopant concentration, a 
small increase in the dopant concentration cause a sharp drop in the resistivity. 
At high dopant concentration, the resistivity approaches that of single-crystal 
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silicon, although it always remains slightly higher (Figure-2). On the other hand, 
carrier mobility in polycrystalline -silicon is lower than carrier mobility in single- 
crystal silicon through out the whole dopant concentration range and is 
characterized by a minimum value at an intermediate dopant concentration [1] to 
[4]- 



N d (cm J ) 


Figure 2: Comparison of Poiysilicon and Single Crystal resistivity 
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Several models have been proposed to explain the resistivity variation with 
dopant concentration in polycrystalline -silicon. In the present work following 
three models have been used: 

(i) The carrier -trapping model 

(ii) The dopant-segregation model 

(iii) The modified model 

Each model explains the electrical properties of polycrystalline -silicon. 

2.2.1 The Carrier-Trapping Model 

This model postulates that the dopant atoms are uniformly distributed through 
out the material and that conductivity is limited by carrier trapping at the grain 
boundaries, where the trapped carriers create potential barriers, which impede 
the transport of free carriers between the grains [1] - [7], [20] and [28], 


Grain 


Grain 


Grain 



Wg b W gb 


Figure 3: A schematic representation of potential barriers 

The high resistivity observed at low dopant concentrations is explained by the 
trapping of most of the carriers at the grain boundaries, leaving few free to 
contribute to the conduction. As the dopant concentration is increased, the 
number of the trapped carriers increases and eventually approaches saturation. 
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Upon further increase in dopant concentration, the number of trapped carriers 
does not increase appreciably; consequently, the potential barriers decrease, 
resulting in a sharp reduction in resistivity [1,4]. Finally, at high dopant 
concentrations, the potential barrier at the grain boundaries become very small 
and no longer limit the conductivity of the samples and the properties of the 
materials approach those of the grains. 

The reason behind carrier trapping is that since the atoms at the grain boundary 
are disordered, there are a large number of defects due to incomplete atomic 
bonding. This results in the formation of trapping state. These trapping states are 
capable of trapping carriers and thereby immobilize them. This reduces the 
number of free carriers available for electrical conduction. After the mobile carrier 
is trapped the trap becomes electrically charged, creating a potential energy 
barrier, which impedes the motion of the carriers from one grain to the 
neighboring grain, there by reducing their mobility. Based on this model, for the 
same amount of doping, the mobility & carrier concentration of a polycrystalline 
semiconductor would be less than that of a single crystalline material. 


2.2.2 Dopant Segregation Model 

In this model grain boundary acts as sink for impurity atoms due to impurity 
segregation at the grain boundary. Consequently the amount of impurity in the 
crystallite is reduced, which leads to a much smaller carrier concentration than 
the uniformly distributed impurity concentration. The carrier concentration does 
not approach that of the doping concentration until the grain boundary is 
saturated with impurity atoms. It is also suggested that segregation of impurity 
caused the grain interior to have higher resistance than the grain boundaries 
[18]. It has been shown that segregation of boron at the grain boundary is 
significant only at extremely heavily doped concentration of silicon. No 
segregation was observed for doping as high as 1.3 x 10 20 cm' 3 [2], [7], [8] & 
[24]. If the reduction 
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of carriers is the result of impurity segregation at the grain boundary, it is 
expected that the carrier concentration reduction would depend on the impurity 
element. It was observed that both boron & phosphorous behaved similarly in 
polysilicon. 

In this model variation in the resistivity is explained by the hypothesis that low 
dopant concentration most of the atoms segregate to the grain boundaries, 
leaving few to contribute to conduction. As the dopant concentration is 
increased, more atoms remain inside the grains, and thus the resistivity 
approaches that of single -crystal silicon. 

Although both the models were successful in explaining the variation of resistivity 
with total dopant concentration, the carrier-trapping model has been more 
generally accepted because the dopant segregation model cannot explain either 
the temperature dependence of the film resistivity or the minimum in Hall mobility 
observed at intermediate dopant concentration [4]. However, neither the trapping 
model nor the segregation model predict any dependence of the electrical 
properties on the high-temperature processing or on the type of dopant species. 


2.2.3 Modified Model 

In this model [4], as the atoms are arranged in a periodic manner forming small 
single crystal, while the grain boundaries are composed of disordered atoms & 
contain large number of defects due to incomplete bonding. The high 
concentration of the defects and dangling bonds at the grain boundaries is the 
origin of the trapping states and the dopant segregation sites around which the 
carrier trapping and dopant segregation model were developed [1] to [9], [20]. In 
these model the grain boundary regions have always been assumed to be very 
narrow compared to the size of the crystallites, thus only the effect of the 
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Segregation sites and trapping states at the grain boundaries on the number of 
active dopant atoms and free carriers have been taken into consideration. 
However, due to highly disordered nature of the grain boundaries and the fact 
that they introduce discontinuities in the otherwise periodic single -crystal nature 
of the crystallites, they should be expected to play an additional direct role in 
determining conduction mechanism in the carriers. 

In order to model the effect of the grain boundaries on the Electrical properties of 
polycrystalline -silicon, some basic assumptions about the properties of the grain 
boundaries have to be made. 

The nature of the grain boundary material can be considered to be between that 
of a completely ordered single crystal and that of highly disordered amorphous 
material. The optical energy gap of amorphous silicon has been measured by 
many researchers [14] -[16], [25] -[27] and was found to be about 1.5 -1.6 eV, 
much larger than the energy gap of single-crystal silicon. Therefore, the energy 
gap of the grain boundary material can be expected to be somewhat larger than 
the energy gap of the single -crystal material within the crystallite. 

The high concentration of defects and dangling bonds at the grain boundaries 
provide ideal sites at which each donor dopant atom at the gain boundary could 
have all of its five bond saturated. Similarly, the three bonds of an acceptor atom 
could be accommodated. Consequently, the dopant atoms at the grain 
boundaries have no weekly bound carriers, and the Fermi level is pinned near 
mid gap at the gain boundaries [14-15] and [28-30]. 

As a result of the above assumption, a hetero-junction is formed at the interface 
between a crystallite and a grain boundary, with the grain boundary material 
behaving as intrinsic wide -band-gap semiconductor. Therefore the effect of the 
grain boundaries can be modeled by potential barriers, shown in Figure -(4). 
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behaving as intrinsic wide -band-gap semiconductor. Therefore the effect of the 
grain boundaries can be modeled by potential barriers, shown in Figure -(4). 


Grain Boundary 

- - ++++-4 

- - +++++ 

- - ++++-+ 

- - ++++-( 

— ++++4 


Depleted Grain 


E c 





X=0 

Figure 4: Energy band diagram for zero applied voltage 


The height of this barrier qfa, relative to the Fermi level, is equal to E g t/2 , where 
E gb is the band gap of the grain boundary material, and the width of each of 
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these potential barriers W gb is approximately equal to the width of the grain 
boundary. Thus in orders to move from one crystallite to other the carriers have 
to either tunnel through the barrier or be sufficiently energetic to be thermally 
emitted over the barrier. 


The potential energy in the depletion region qV (x) relative to the valence band 
edge E v in the neutral region of the grain is given as a function of distance x by 
reference [1], [2], [4] and [20]. (Appendix A) 
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2.3 W.K.B Method 

The grain boundary material has a resistivity much higher than that of the grains 
due to its intrinsic and highly disordered nature. Thus conduction mainly takes 
place by carrier motion from one grain to another by tunneling though and/or by 
thermionic emission over the barriers at grain boundary. 

The one-dimensional time -dependent tunneling probability D (E x ) for a barrier of 
height qV (x) is given in the W.K.B approximation using references [1 1-12] 
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D(E X ) = exp 





Where E x is the energy component of the incident carrier in the X-direction with 
zero energy taken at the top of the valence band in the neutral region of the 
grains, and , the effective mass of the carriers is assumed to be energy 

independent. (Appendix-B). The validity of the W.K.B tunneling probability 
approximation has been studied theoretically and experimentally by [17,39], and 
was found that. Polysilicon satisfies all the conditions (small grain boundary 
width, depletion potential and grain boundary potential) shown in [39] to use 
W.K.B technique. 


2.4 Fermi-Dirac Statistics 


Fermi-Dirac statistics give the distribution of the holes having energy E by the 
following equation. 

N(E) = [l-f(E)fi 
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N(E) = w 
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1 + exp 
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Where E F is the Fermi energy & w is the weight factor. Above equation show the 
number of holes [13]. 
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Chapter 3 

Analytical Calculations 


3.1 Hole Density Incident On The Barrier 

For the calculations of the current density (J) the hole density incident on the 
barrier per unit area & per unit time is required. To calculate hole density Fermi- 
Dirac statistics is used, according to which, out of N holes N (E r ) are in the state 

E r and weight w is given by 


Or 


f 
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N(E) = w 


1 - 


" {E-E f ) 
1 + exp - — — 


kT 


jy 


N(E) = w 


VL 


1 + exp 




kT 


In the phase space-weighting factor is given by, [14]& [19]. 




dp x dp y dp 2 


(la) 


(1b) 


(2) 


Here V is the volume in spatial coordinates i.e. V = dxdydz 
From equation (1) and (2), 
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2Vdp r dp^dp. 

N(E X , y ,;)= - 3 ^ - ' 


1 + exp 


(Ef-E) 


Number of holes moving in x-direction in range dp x per unit area and time, 


9 j 00 00 

N(E X )= \dp y J 

hm k -00 -00 


£ -E 

1 + exp — - 

kT 


The total kinetic energy E r in the x-direction is given by 


E r ~ ^r(Px + P 2 y + Pi) 


Solving these equation using polar co-ordinates in Y-Z plane 


••• pI + pl = pi 


dp y dp z = p dp dO 


So the equation become (4) 


2 p dp 2n<x> 

v( E x ) = -j-r- J J- 
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And its solution is given by, 


N(E r ) 


A; -~e f 

Ait J 2m h 

kTp v dp v ln^ 1+exp — — 
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Converting the equation in terms of energy and taking the hole kinetic energy 
w.r.t valence band. 



Where %=E F -Ev. 

This is the final expression of the hole density entering at the potential barrier per 
unit time, per unit area (Appendix-C) 

3.2 Current Density (J) Produced By Charge Carriers 

Let V 0 be the voltage drop across the complete grain boundary region. This will 
be drop in voltage across the grain -boundary. Barrier and the two depletion 
regions on either side of the grain boundary, are as shown in the Figure -(6) such 
that. 



Figure 6: Energy band diagram near a grain boundary at an applied voltage 
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Now the number of electrons tunneling from left to right Nlr and number of 
electron tunneling from right to left Nrl can be calculated by integrating 
D(E X )N(4,E X ,T) over all energies & using (2.3) with equation (8). 


4ft * °? . 
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( 10 ) 


In general D L r (Ex) and D RL (E x ) are function of different voltages dropped 
across the grain boundary regions on its sides so that D L r(Ex) ^Drl(Ex). 

Thus the net current density (J) through the grain boundary is given by 


J - q [Nlr~NrJ 


So putting the Nlr and Nrl from equation (10) we get, 
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Now to solve equation (11), we need tunneling probabilities D L r(E x ) and D RL (Ex), 
both of which depend on barrier qV (x). From the Figure it is clear that q V (x) can 
have different forms depending on the energy of the carriers. 

So we have three different conditions. 

(i) E x <qV (x) (TFE): For the carriers, which have energies less than the 
energy required to over come the potential barrier in the depleted regions of the 
grains, the potentials barrier qV (x) is defined as follows, 

qV (x) = 0 for\x\ >W + W gb /2 

qV(x)= +W a /2)-\x^ forW lt IK\x\<W + W sl /2 

qV(x) = qj HA< W *I 2 

In this case thermionic field emission is possible (TFE) 

(ii) qV b <E x <q0 (TFES): If the carriers have sufficient energy to over come the 
parabolic barriers in the depleted regions but not sufficient energy to over come 
the grain boundary barrier qfa then the potential barrier qV(x) through which the 
carriers have to tunnel is given by 

qV(x) = 0 for \x\>W gb /2 
qV(x) = q# for |*| < W gb j2 
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In this case Thermionic field emission with scattering (TFES) from the grain 
boundary is possible. 

(iii) E x > q<f> (TE): If the carriers are sufficiently energetic to over come the 
grain boundary barrier q^, then qV(x) = 0, D(E X ) = 1, and current transport is by 
the Thermionic emission over these barriers. So in this case thermionic emission 
(TE) is possible. 

In general all the above conduction mechanisms occur in parallel. However, at a 
given temperature or with in a small range of temperature, carrier transport is 
dominated by one or two of these process. The present work is restricted to a 
practical range of temperatures and the following limitations have also been 
imposed while obtaining an analytical solution. 

[i] The solution of current density is restricted to temperature at which the 
carriers have sufficient energy to over come the depletion region 
barriers tunnel through the much thinner grain- boundary barriers 
becomes more significant. 

[ii] The average applied voltage per gain V 0 is much less than the height 
of the grain boundary barrier ( V 0 « q<j > ). 

[iii] The maximum barrier height in the depletion region of the grains is 
much less than the height of the grain -boundary barrier ( Vb< <q<P) 

[iv] Complete ionization of the dopant atoms in the grains is assumed so 
that the carrier concentration in the neutral regions of the grains is 
equal to the dopant concentration there. 

[v] The single -crystal energy band is applicable inside the crystallites. 

[vi] One-dimensional majority -carrier transport is assumed, any 
contribution from the injection of minority carrier is neglected. The 
equivalent effective mass of single crystal silicon is used in the barrier. 

At intermediate temperature only a few carrier have energies higher than the 
maximum height of the grain boundary barriers. Thus the contribution to current 
by Thermionic emission of carriers over the grain boundary barriers can be 
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neglected. Also the grain -boundary barrier are much thinner than the barriers in 
the depleted regions of the grains. Thus the probability of tunneling through the 
grain -boundary is higher than passing through depletion region & grain 
boundary barrier. Thus conduction is mainly controlled by carriers tunneling 
through the grain boundary barriers after being thermally emitted over the 
depletion region barriers; thus carriers whose value of E x are in the range qVb < 
E x < q</>. By including the above approximations equation (11) can be written in 
the form. 


j = 




The above equation can further be simplified by considering only the carrier 
distributions ln(1+X)&X. so the equation become, 


■k 


AK m h kT 


ft<t>+qV g bj1 
\)V b +qV g b/2 ^ ^ 

Wgbl 1 j "I 

< exp 

Vb-qVgbl 2 


U 

r 


ExH+^+qVR 


kT 


■E>LR(E x )dE x 


I 

} 


DRL i^x) 


(15) 


As current flow is predominantly due to carriers whose value of E x are only 
slightly greater than the maximum height of the barriers in the depletion regions, 
thus the tunneling probability D(E X ) can be simplified by expanding it using Taylor 
series around qVb .( Appendix-D) 
s x = qV b - E x 


.: In D(E x ) = ~[bi+ cie x + fis x + 


■] 


(16) 



Analytical Calculations 


21 


Where, 

b i = (17) 

c ] =!^-)[2m;( q V(x)-qV b )pdxj 

For applied voltage much less than height of the grain boundary barrier 
q</> (assumption [ii] made earlier), bi , Ci and fi are independent of the applied 
voltage so that 

Dlr(E x ) = Drl(E x ) = D(E X ) (18) 

The constants bi and Ci can be calculated by solving (17) using the appropriate 
barrier qV(x). For carrier in the energy range qV b <E x <q</>, the values of b i, Ci are 
given by (Appendix-E). 


6,= T w *>k m 'kT { ai-i v i>f 12 

c\=^- w gb(?>”hf 2 (<it-<if r b)~' /2 <19) 

J 

The quadratic & higher order terms in the expansion of D(E X ) can be neglected 
using assumption [iii] i.e. qV b « q<f> ,as they will contain the terms of 

{q$-qV b )T ,(q<fi-qV b )~2 etc. which become less & less as the negative power 
increases. 

So the value of the D(E X ) can be taken as follows; 

D(E X ) = exp[-(bi + cis x )] (20) 


Using equation (20) in (15) the final expression of the current density given by; 
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j 


= <7 




^v[-{b\+c\E x )]dE x 


>ex.p[-{bi+c\s x )]dE x 


(21) 


After solving the integration (appendix-F) we get, 



(22) 



q v L + 
e' 




(23) 


Since according to our assumption q<p » qV b (for TFES)), thus ignoring terms of 
q<j> in the exponential we get, 



( 24 ) 
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Equation (24) shows that current density is a function of both the voltage drop 
across the grain boundary barrier and the voltage drops across the depleted 
regions on either side of grain boundary barrier. For the tunneling to be occurring 
the applied voltage V 0 should be much less than the height of the barrier in the 
depletion region (Vo « V b ) under this condition voltage drop on left and right side 
is almost similar [4]. 

• • v L =v R = v d 


Vdis half of the voltage drop across the depletion region of the grains. 
So the new from of (24) is, 


J = 


0-q«') 


4?r *, 2 t 2 

—mhk T 


) 


~4/kT -by - 
e e 


■qVb/kT 


2 Sink 


qVd+cykTq Vgb /^ 


kT 


(25) 


Now the voltage drop across the grain boundary (V gb ) is given by (Appendix G) 
i.e., 


gb 


WgbNp 

Qt 


2V, 


Also 


Vo = V gb+ 2V d 

W gb N G 


V n = 


Qt 


2V d+ 2V d 


V 0 =2V d 


Qt 


+i 


Using equation (26a) & (26b) in (25) we get; 


(26 a) 


(26 b) 


J = 


2q ( 4 * * 2 l\-^ kT qV ^ kT 


e e 


q y„ 

1 1 NGWgb 2 kT 

Qt 


(27) 
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This is the expression of the current density; It can be seen that the current 
density is dependent on the applied voltage, trapping density & barrier height of 
the grain region. Assuming that applied voltage V 0 across the composite grain 
boundary is much less than AT, the assumption reduces equation (27) to 


J 


2q 

~ (l-qAT) 



/ \ 

;\~ q¥blkT 

\}P / 

1 qVo 

S' 

to* 


(28) 


Now as we know that, 


J — Ggb s gb ~ gb 



Therefore the conductivity of the composite grain boundary region a 9 b is given by 



f \ 

( 4* « 1 Sl kT -*l -* V bl kT 

\ h 3 mkkT r e e 

l+c\kTNQ Wgb /Q t 
,/GWgb ' 

l Qt ) 


(29) 


where L gb =2 W + W gb is the length of the composite grain boundary region. 

The average resistivity can be obtained as follows; in general the voltage drop 
across a single crystal V a is divided between the composite grain-boundary 
region (the grain boundary and the adjacent region in the grains) and the grain 
itself. 


V^Vo+Vg 

Where V g is the voltage drop across the undepleted grain. 


P = 

P = 


Va 

JL 


s 


v„ _ V* 


— * O j v ... 

~ JL JL~ UL, 


8 






(30) 
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3.3 Calculation of Resistivity & Mobility for Various Doping Conditions 

Current density or the conductivity expression was derived under the assumption 
listed before but with out any limitation on the dopant concentration. However, 
the width of depleted regions in the grains due to carrier trapping is a function of 
the dopant concentration. Therefore, the dopant concentration indirectly affects 
the conditions under which the above theory is applicable. 


3.3.1 At High Dopant Concentration 

From equation (29) we have 



Assuming the complete ionization of the dopant atoms within the grain the carrier 
concentration in the grain is given by. 


n = N v exp 




z£ 

\kT j 


^ 2mn h kT^ /2 

C * j 


exp 


r z£] 

\kT J 


N r . 


With this condition equation (29) becomes as follows 


q 2 Lgb N G 

r gb = (l-qW) 


1 




> 


/ 2 qV VkT 


Qt 


NqWgb 


+c\kT 


u — 


Qt 


N G W gb 


(31) 


(32a) 
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And, 


q 2 NQ 

J TFES - Ti=q Wl 


( , V 

f \ 

® +c\kT 

N G W gb 

\ * 1 e e 
'iTm^kT ' 

i + a ' 


N G W gb J 


(32b) 


Equation (32a) is the conductivity equation for highly B-doped polysilicon, which 
depends on grain boundary width and depletion potential barrier. It is a 
transcendental equation for W gb ; in the next chapter calculation of W gb is done 

numerically using some experimental values and fitting parameter. So starting 
with equation (19), values of constants b-i & c? depend on three parameter 
W gb ,<7{0 and qV b . In the present work we have calculated the W gb using the 

experimental values of q<f>an<\qV b so for a particular sample q<j) and qV b will be 
constant. Thus bi and ci are taken as function of W gb as follows. 

b-t-bWgb & Ci~cWg b (33a) 


Where 


Also define 



Using equation (33a) in equation (33b) we get, 



(33b) 


( 34 ) 
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3.3.2 Low Dopant Concentration 

At low dopant concentration the height of the depletion region barrier qVb has a 
maximum value approximately Eg/2, where E g is the band gap of the single 
crystal silicon (l.leV). Height of the grain boundary barrier q<f> g b- E g t/2 is 
expected to be of the order of 0.6 —0.8 eV (Section IV, i). Thus at low dopant 
concentration, qVb may be comparable to q<p. 

So for low dopant concentration q<j> &qVb thus the assumption [iii] has been 
violated and hence [i]. It must be noted that assumptions [iii] & [i] have been 
used in many places to derive equation (29). Since depletion barrier height is 
comparable to grain boundary barrier, the tunneling probability becomes one & 
thermionic emission (TE) will be dominant over thermionic emission field 
emission with scattering (TFES). Therefore, at low dopant concentration, 
conduction properties can be calculated by thermionic emission (TE) and we 
cannot calculate the properties using thermionic field emission with scattering 
(TFES). 


3.3.3 Thermionic Emission (TE) 

As mentioned above, at low doping concentration the depletion barrier is 
comparable to potential barrier across the grain boundary. Hence tunneling 
probability is unity. Therefore thermionic emission becomes dominant. To 
calculate current density due to TE we simply put 1 for tunneling probability in 
equation (15). Hence 


D(E X ) - exp[-(bi + cis x )]=1 orbi,ci- 0 ( 35 ) 

Hence in equation (25) given below, 


j = 


t 


(l-q«') V 


^mlkh 2 


) 


-4/kT -by -qV b /kT 

e e 2 Sink 


qV d +c\kTq 




kT 


(25) 
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After putting bi = 0 & Ci= 0 equation (25) becomes, 


[ 4 n * 2 2 \S /kT ~ qV b! kT qV d 

J = q\-^m h k l T l fe e e 2Sinh SIS. 

\h J ' kT 


v V Q «kT Consequently V d «kT and assuming the complete ionization of the 
dopant atoms within the grain the carrier concentration in the grain is given by. 

xr (-£' j2am* h kT\ ( -%') sr 

« = exp =2 — exp — « N G 

\kT J I h I y 


With these values equation (36) gives the thermionic current density as, 


J TJ? ^ N Ct * 

TE U \27m h kT 


-qV b /kT 


Using the relation of equation (26b) 




-qV b /kT 




The conductivity is given by, 


Hence 


q 2 LN G i_ 

W gb NQ J [ 2nm b kT 


-qV b /kT 


(38a) 


Vg f 

Qt 

q 2 LN G 


* 'll/ 2 qV b/ kT 

2nmfokT ) e 


(38b) 
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And 


qL 


TE 


W gb N G 


Qt 


+1 


1 


iTmpJcT 


1/2 


-qVb/kT 


(38c) 


This is the conductivity equation due to thermionic emission. We can further 


simplify this equation, by assigning 


W eb N c 


Q, 


«1 


••• °TE = 9 LN G 


27mhkT 


1/2 


-qV b jkT 


(38d) 


This equation is same as obtained by Seto [2], Lu [7], Baccarani [3], though they 
neglect the grain boundary effect completely, equation (38d) better explain the 
experimental results at low doping concentration. 

Equation (38c) is the mobility equation at low doping concentration. We will see 
in chapter 4 that it will show the deviation from the experimental values at high 
dopant concentration. In the following derivation we have shown that polysilicon 
mobility reaches to single crystal range at high dopant concentration. We have 
equation (37), 


J TE q N G 


xl/2 


1 


Irnn^kT 


-qVblkT v 0 


Qt 


Say G = 


'W&N 0+ ^ 


V 


Qt 


j 


••• J T£ = <7 M G 


1 


Iron^kT 


1/2 


-qVb/kT 
e V G 


(37a) 
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Current density due to undepleted region. 

J q— o’qEq = qN G fj. y (37b) 

mono 7 — ~ — r 
(L-2PV) 

where n mono is the mobility of single crystal. 

Current density must be maintained at any point of the sample, therefore 
equating equations (37a) and (37b). 

J G ~ J TE 


qN Q M 


V 


g 

mono-, — 2 — r= q N - 
(. L-2W ) H G 


1/2 


-qVb/kT 
e V G 


IrnnfjkT J u 

If V a is the applied voltage V a =V 0 +V g , Eliminating V 0 we get, 

qG 


V s = 


{inmlkT )' 


l '/2 


-M/kTy 


Mmono qG _ e -qV b /kT 


(L~2W) ( 2m l k Tf 


Final expression of current density from equation (37b) is given by using this 
potential across the grain (Vg), 


qN c 


V* 


qG 


J G = 


(L- 


-Rh/kTy 


qG 


(L-2W) ( 2 mnlkT) 


,-qyjkT 
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Dividing both sides by L and rearranging the equation. 


(&a) 

(K) 


'{L-2W) | (2 nm[kTf MtT ~ 

^ G ) 

UJ 


_ M mano L qLG 


If Eext is the external applied field and V d is the drift velocity of a hole. Using 
equation (38c) we get, 


( e 1 

^ ext 



, 

+ 

^ , 
1 
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_ t^mono L f^TE J 
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Ae# 


1 (L-2W) [ 1 

Mmono ^ Arc 


At high dopant concentration 2W« L 


1 


Ae# 


1 1 
1 

_ mono f^TE _ 


Or fi eff = 


Mmono f^TE 

(A m0 „ 0 +Are) 


(39) 


This equation show how effective mobility related with singe crystal and 
thermionic mobility.. In chapter 4 plot is shown which show the difference of 
A ^ , Arc an d the experimental values at high dopant concentration. 


3.3.4 Combined Effect of TFES & TE 

While deriving thermionic emission we have neglected the effect of scattering or 
tunneling through the barrier, though we cannot neglect the tunneling effect 
completely there is always some probability for the tunneling so their combined 
mechanism should be considered. Considering both TFES and TE effects the 
current density from equations (32b) and (37) is given by, 
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J Total J TFES + ^ TE 


J = g N 


G 


\ 1/2 


V (TFES + TE) 


(40) 

(40a) 


^ 27tmfokT ' 

Now the final expression of the resistivity is given by equation (30) 


P L P *+L P » 


^ gb 

p-~Y"Pgb + p g 


2 W W gb 


(30) 


Using the relation J = 


P = 


PgbLgb 


and equation (40a) in equation (30) we get, 


,V2 


1 


ilrnnlkl) 

q 1 LN „ \TFES + TE) 
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g 


1- 


2 W W gb } 


(41) 


This is the average resistivity equation of the polysilicon. Using this equation 
calculation of resistivity at low and high dopant concentration is done. 


3.3.4. 1 At High Dopant Concentration 

At high dopant concentration only TFES is effective and the contribution by 
thermionic emission is negligible, so average resistivity is given by thermionic 
field emission with scattering (TFES) only. So putting equation (34) in (30) we get 
the average resistivity at high doping concentration. 


_ {^Tcm\kT )' 


l '/2 


(\-cW tb kA\+— Q— 
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(42a) 
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If we also consider the TE with TFES the equation changes to, 


( 2 jrmlkT )' /2 qVb/kT 
q 2 LN a 



(l-cW gb kT 1 1 + 


Wo) 


vr^» +1 ' 


V 


Q, 
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( x 2 w 
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And 


1 

M = WgP 


(42b) 

(42c) 


This is the resistivity equation at high dopant concentration, which depends on 
grain boundary width, grain size, trapping density, potential barrier, temperature 
and resistivity of the crystallites. Equation (42b) is used for the numerical 
calculation of grain boundary width using experimental results in the next 
chapter. 


3.3.4.2 At Low Dopant Concentration 

Contribution to the current density by both TFES and TE at low doping 
concentration is the sum of both the currents in the grain 

J Total = JtFED + JtE 

As explained in section 3.3.2 at low doping concentration thermionic emission is 
the dominant phenomenon, so current density due to thermionic field emission 
with scattering ( J TFED ) can be neglected. 
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J Total J TE 
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(43a) 

(43b) 


(43c) 


Using Equation (43b) effective resistivity can be calculated at low dopant 
concentration. 


P = 


q 2 N G L 


(43d) 


In the present work special interest is shown to calculate average carrier 
concentration at low doping concentration, to calculate average carrier 
concentration we have used the following equation. 


No 


1 

qpfi 


(43e) 


These equations (43a-d) are used to calculate current density, mobility & 
resistivity for various low doping conditions and the results are compared with 
experimental values in the next chapter. 




Chapter 4 

Numerical Calculations And Results 


4.1.1 Numerical Calculation of W gb 

In chapter 3 derivation of resistivity is done analytically (equation 42). In the 
present section numerical calculation of grain boundary width is done. The 
experimental data have been used to numerically calculate W gb . To solve 
equation (42) numerically we need the depletion width, which can be calculated 
using the following formula for single crystal silicon [10] using the diode 
approximation between depletion width and grain boundary (assumption [v]). 


W = 


2 sV, 


qN c 


1 + ^ 


V 


N, 


' J 


(44a) 


where N c &N t are the charge densities within grain and grain boundary. Charge 
density inside the grain boundary is given by carriers trapped inside the grain 
boundary, i.e. N, 

Now the two cases arise for the depletion width. 

(i) Depletion width at very high dopant concentration 

(ii) Depletion width at moderate dopant concentration 


In case (i), at high dopant concentration depletion barrier height is extremely low. 
Under this case V b = 0 and consequently W«0 from equation (44a). With these 
values in equation (42a) gets modified as follows. 
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p = 


( Imri^kT y 2 
q*LN a 


{\-cW gb kT 

^ G w gb 



(45) 


Equation of average resistivity only due to thermionic field emission with 
scattering from the barrier (TFES) is taken for the grain boundary width 
calculation (W g b) as at high dopant concentration TFES is a dominant 
phenomenon, so we had neglected contribution due to TE. To solve equation 
(45) we need the following experimental data. 


S.No. 

Data Taken from Experiments 


Physical parameter 

Value 

1. 

Q, 

3.34xl0 12 cm~ 2 

2. 

N g 

2.6x10 19 cot* 3 

3. 

P 

4.2xl0 -3 Q.cm 

4. 

Pg 

5.18x10~ 3 Q.<7m 

5. 

q<l> 

MeV 

6. 

L 

200 A 0 

7. 

T 

27 °C 

8. 

£ 

1. 0433x1 0’ 1U 

9. 

ml 

.23m 0 =2.093 xl0“ 3 ’ 

10. 


480c/n 2 F _1 sec' 1 

11. 

K 

i 

1.38x10 19 kg.cm 2 .sce~ 2 K~ 1 


Table 1: The Experimental data from [2,4 & 7] for very high doping 


Calculated value of b & c using equations (33b) and values given in the Table 
above 
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1. 

b 

3.9910xa10'cm" 

2. 

c 

1.8897xl0 22 %'W 3 sec 2 


Table 2: Constants calculated 


Using the tabulated values, transcendental equation (45) can be solved 
numerically. The Newton-Raphson iteration method (Appendix-H) in the Mat Lab 
application software was used [34] to calculate W gb . W gb of ~ 11.86A 0 was 
obtained for high doping. 

Now in case (ii) when the doping is moderate, barrier height and depletion width 
are not negligible, so these have also been calculated. 
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2 eV b 
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1 + ^2- 
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2 sV b 

qN G 

Lwgb) 
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(44b) 


and equation (42a) 

(l-cW ^kT )| 


_ ( 2nm\kT )' 


,1/2 
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Q, 


W .N 


c J 


q LN 


Q, 


K N c w * 


-+ ckTW 


gw# + ^ 2 IT W 0 j 


Po 

(42a) 


Equation (42a) is a transcendental equation, using equation (44b) and following 
the similar procedure mentioned above with the following experimental values 
given in the Table 3 and 4. 
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S.No. 

Data Taken from Experiments 


Physical parameter 

Value 

1. 

Qt 

3.34x10 12 c/w -2 


N g 

5x10 1! W 3 

3. 

qV b 

0.0335eF 

4. 

P 

0.85fi.cw 

5. 

Pg 

0.02Q.cjm 

6. 

q<t> 

MeV 

7. 

L 

200 A 0 

8. 

T 

27 °C 

9. 

6 

1. 0433x1 0' 1U 

10. 

ml 

.23m 0 = 2.093 xlO" 31 

11. 


60 cm 2 V~ l sec -1 


K 

1.38 x IQ 19 kg.cm 2 .sce~ 2 


Table 3: The Experimental data from [2,4& 7] for very high doping 


Calculated value of b, c using equations (33b) and values given in the Table 
above. 


1. 

B 

3.8884xl0W l 

2. 

C 

1 .9395 x 1 0 22 kg~ l cm~ 3 sce 2 


Table 4: Constants calculated 


Using these values in equation (42a) W gb comes out to be &11.89A 0 . It can be 
seen that change in the value of grain boundary width is negligible with the 
change in dopant concentration. So grain boundary width is not a function of 
dopant concentration. Similar calculation was done for grain size L=1220A° and 
the trapping density Qt= 1.9x10 12 cm' 2 at room temperature the grain boundary 
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width comes out to be 27.35A 0 for very high doping and 27.87A 0 for moderately 
high doping concentration. 


4.1.2 Resistivity at High Doping Concentration 

Resistivity at high doping concentration is governed by equation (42), 


p = 


( 2nm\kT )' 
q 1 LN\ 




‘iVbpT 
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{N G W gi 


- + ckTW 
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~bW 
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(l -cw gb kT ) 
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Q, 




'* w * N cj 


Q, 


+ 1 
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P{tfes+te) 




Pg 


If we put the numerical values in the equation, we find that at high doping 
concentration, resistivity is of the range of single crystal as shown in Table 5, and 
the first term in the above equation makes negligible contribution to resistivity. 
Thus at high doping concentration resistivity of grain becomes dominant. 



L=200A° 

S.No. 

Doping 

Concentration 

(cm- 3 ) 

Experimental 

Resistivity 

Resistivity of 
Composite G.B 

Region 

P(TFES +TE ) 

[Q.c m\ 

Resistivity of non 

depleted Grain 

Pg 

[Q.cw] 

1. 

5x10'° 

0.85 

7.4829x1 O' 4 

.846 

2. 


.00278 


.00238 

3. 

5x10 ia 

4.2x1 0" 3 

3.35x10'“ 



Table 5: Comparisons Between Resistivities 
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4.2 Theoretically Calculated Values & Various Plots 

Using equations derived in the previous chapter along with the parameters of the 
Tables 1-4 the different electrical parameters average carrier concentration, 
barrier height, mobility, resistivity and current density are calculated and 
presented in the respective Tables at room temperature. In the present work our 
special interest is to calculate average carrier concentration at low doping 
concentration. We have used equation (43e) to calculate average carrier 
concentration. Therefore we need mobility and resistivity, hence barrier height. A 
brief procedure how these parameters are calculated is shown below. 

(i) Barrier height given in Table 6a is calculated using equations from 
section 2.2.3. Values of Qt, L, q and s are given in Table 1. Trapping 
density 1.9x10 12 cm' 2 is taken for £.=1220A°from[7,20]. After calculating 
barrier height for different doping concentration and grain size, 
Average mobility shown in Table 6a was calculated using the empirical 
equation (43f) from section 4.2.1 and using the data of Table 1. 
Different mobility /u mono ,v TE and /u eff shown in Table 6b. fi mono is taken 

from [21] and Are./V are calculated using equations (38c) and (39) 
respectively. 

(ii) Table (7-8) shows the mobility variation considering contribution from 
TE, TFES, and TFES+TE along with experimental values at two 
different grain size L=200 A 0 and £=1220A° respectively. Equations 
(38c & 39), (42c & 42a) and (42c & 42b) are used for mobility 
calculation due to TE, TFES and TFES+TE respectively. Experimental 
values of Mobilities for the comparison are taken from [2] and 
[20,21 ,24] for £=200 A 0 and L= 1220 A 0 respectively. 

(iii) Table 9 shows the mobility values as a function of 1 1kT. All values 
were calculated using equation (43f) from section 4.2.1 for grain size 
1220A° and barrier potential 0.078eV[20]. 

(iv) Average carrier concentration shown in Table 10 is calculated using 
equation (43e). Resistivity values are taken from [2,7,20] for i.=200A° 
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and from [7,20] for L=1220A° the calculated average mobility values 
shown in Table 6a are used for the calculation. 

(v) Table (11-13) shows the resistivity variation considering contribution 
from TE, TFES and TFES+TE along with experimental values at two 
different grains size L=200A° and L= 1220 A° respectively. Equations 
(38b), (42a) and (42b) are used for the resistivity calculation due to TE, 
TFES and TFES+TE respectively. Experimental values of resistivity 
are taken from [2,7,20] for L=200A° and from [7,20] for L=1 220/4°. 

(vi) Table 14 shows the resistivity values at different grain size. Values are 
calculated using equation (43d) at doping concentration A/=1x10 17 cm' 3 
and trapping density 4.4x1 0 12 cm' 2 [32], 

(vii) Table 15 shows the current density at different operating voltages. 
These values were calculated using equation (43a), other experimental 
data needed are taken from [2], [22]. 


4.2.1 Plots and Corresponding Data of “Barrier Height Vs Doping 
Concentration” and “Mobility Vs Doping Concentration” 


Plot of mobility Vs doping concentration is drawn using equation (43c), to match 
calculated values with experimental results two additional fitting parameters f and 
n are introduced here. These fitting parameters cannot be attributed to 
theoretical inaccuracy; it is due to two possible reasons. One is that the effective 
Richardson constant (effective hole mass) is much smaller than in single crystal, 
which may occur in such a disordered structures as polysilicon. The second is 
that some transmission probability exists when carriers pass through the 
complicated grain boundary by either scattering or recombination. 

So equation (43c) with these fitting parameters is given by, 


a = fm 


Inm^kT 


1/2 


-qVb/nkT 

e 


(43f) 
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f is a scale factor which modifies the effective Richardson, S constant of a single 
crystal when used in such a highly disordered structure. Now equation (43f) 
becomes an empirical equation analogous to resistivity equation in [7,20,21]. In 
Table 6, mobility calculated from empirical equation (43f) is shown and 
compared with the experimental mobility values [7,20] in the following sections. 
Values of the fitting parameters from [7] are given below, 

£=0.12 , n= 6.49 for L=200A° 
f= 0.060, n= 1 .22 for L=1220A° 
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Grain Size=1220 A 


0.0003 


07 


0.0014 


0.002 


0.0029 


0.0071 


0.0143 


0.0214 


Grain Size=200 A 



157.3008 


158.7290 


151.7153 


149.6323 


145.0114 


126.6255 


101.0150 


80.5844 


64.2859 


5.0x10 


7.5x10 


3.83x10 


0.0006 


0.0008 


0.0031 


50.8984 


46.2508 


11. 3.1x10 1B 

0.0885 

9.6334 

■KE333EM 

0.0090 

■■rearmii 

12. 3.2x1 0 1b 

0.0913 

8.8009 

wmmm 

0.0367 

■MMHTTM 

IHHEnifli! 

0.0999 

6.7106 

BESEEM 

0.0699 

1 3.4881 1 


0.1215 


0.1281 



28 


285 


0.0214 


0.4749 


0.3690 


.7592 


13.4257 


21.0900 


26.4330 


20. 

I.OxlO 1 ' 

0.2853 

0.0188 

MEiSODtMil 

0.0086 

39.6892 

21. 

■H333EMH 

0.4423 

0.0001 

MEE5QMH 

0.0043 

45.4480 

22. 

IMB333EMH 

0.2403 

0.0247 

wKEsmm 

0.0029 

47.5476 



0.1267 


0.0692 


0.0149 


0.0014 


1.9793 


66.0428 


102.3861 


151.9199 


155.2872 




48.6335 




Table 6 a:Doping Concentration, Barrier Height and Mobility 


Doping Concentration 
(cm' 3 ) 



(cm /V.sec)[21] 


562.34 


554.30 


453.15 


392.41 




36 


294.27 


165.48 

17.720 

16.0099 

133.35 

66.040 

44.1685 

100.00 

102.38 

50.5895 


93.05 


86.59 


64.93 


6.23 


54.63 


53.85 



155.28 


39.9887 


Table 6 b: Doping Concentration, Barrier Height and Mobility 
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4.2.2 Plots and Corresponding Data of “Mobility Vs Doping Concentration” 
For TFES+TE, TFES, TE and Experimental Values (200 A°Grain Size) 




Dopant concentration (cm' 3 ) 

Figure 9b: Mobility Vs. Concentration for TE+TFES, TFES, and Experimental Value, L=200A° 
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S.No 

Dopant 

Concentration 

[cm -3 ] 

For Grain Size= 200^4° 

Experimental 
Mobility 
\cm~ 2 V~ l sec -1 ] 
[2] 

TFES+TE 
Mobility 
[cm -2 F -1 sec -1 ] 

TFES 
Mobility 
[< cm~ 2 V~ x sec -1 ] 

TE 

Mobility 
[ cm~ 2 V~ x sec -1 ] 


1.00X10 lB 




45.70882 


— 

52.48075 



39.81072 


1.20X10 

41.68694 



■ESSRSM 


ikbseqmi 

16.2181 



12.8825 


IHSiMl 

5.4881 





hh 




0.4749 


■HH 


1.36901 

0.36901 

1.148154 


2.00X1 0 18 

1.230269 

2.7592 

1.023293 

2.844461 




5.623413 

2.630268 

19.6336 


7.50X1 0 18 


15.13561 

9.120108 

31.62278 


1.00X10 19 1 

15.4257 

20.41738 

13.4257 

39.81072 


2.50X1 0 19 

23.09 

31.62278 

21.0901 

56.23413 


5.00X1 0 1a 


39.81072 

26.433 

63.09573 


7.50X10 19 1 


52.48075 

35.48134 

66.06934 

■ 


41.6892 

56.23413 

39.6892 

69.1831 


Table 7: Concentration and Mobilities due to different effects 
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4.2.3 Plots and Corresponding Data of “Mobility Vs Doping Concentration 
For TFES+TE, TFES, TE and Experimental Values (1220A°Grain Size). 


Dopant concentration (cm ) 

Figure 1 0a: Mobility Vs Concentration for TE+TFES, TFES, and Expt. Values for L=1220A 


0 TFES 
* TFES+TE 
Exi 


Dopant concentration (cm' 3 ) 

Figure 10b: Mobility Vs Concentration for TE+TFES, TFES, and Expt. Values for L=1220A 
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Dopant 

Concentration 

[cm- 1 ] 

For Grain Size=1220A° 



TFES 

Mobility 
\cm~ 2 V~ x sec -1 ] 


■ 

■ntrinuM 

64.28590 



66.2859 


■ HMiniHi 

34.67369 



28.0361 


3.50X10 10 

10.54090 



HM 


nm mm 

2.137962 



1.737801 


9.00X1 0 1D 




0.645654 

MM 

KHQM 




0.316228 

|| 

i^MM 


0.000278 

0.0001 

0.000278 

HIM 

BgESEMM 

0.123027 

0.097724 

0.380189 

0.097724 


4.57X1 0 1 ' 

2.630268 

1.819701 

5.888437 

2.570396 


EEHM 

40.73803 

56.23413 

36.30781 

56.23413 


eesemm 

72.4436 

117.4898 

58.88437 

112.2018 

|| 

mnnMH 

91.20108 

151.3561 

75.85776 

138.0384 

mm 

EEMMMi 

91.20108 

144.544 

75.85776 

141.2538 


Table 8: Concentration and Mobilities due to different effects 
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Average Carrier Concentration (cm‘ 3 ) Average Carrier Concentration (cm 3 ) 


4.2.5 Plots and Corresponding Data of “Doping Concentration Vs Average 
Carrier Concentration” 



Calculated Doping Concentration 



Figure 12b: Average Carrier Concentration Vs. Doping Concentration at Low Dopin a le 
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Table 10: Doping Concentration and Average Carrier Concentration 
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4.2.6 Plots and Corresponding Data of “Resistivity vs. Doping 
Concentration” For TFES+TE, TFES, TE and Experimental Values 
(200A°Grain size). 



Figure 13a: Resistivity Vs Doping concentration for TE and Experimental Values for L=200A° 



Figure 13b: Resistivity Vs Doping concentration for TFES+TE, TFES and Experimental Values 
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S.No 

Dopant 

Concentration 

[cm -3 ] 

For Grain Size= 20(M° 

Experimental 

Resistivity 

[p.cm][2,7] 

TFES+TE 

Resistivity 

[Q.cw] 

TFES 

Resistivity 

[Q.cm] 

TE 

Resistivity 

[Q.c/m] 

1. 
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2. 
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9. 

niiiimif 




HriUcVTCiisaa 


T able 1 1 : Concentration and Resistivity due to different effects. 
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4.2.7 Plots and Corresponding Data of “Resistivity vs. 
Concentration” For TFES+TE, TFES, TE and Experimental 
(1220A°Grain size) 



Doping 

Values 


Dopant concentration (cm' 3 ) 

Figure 14a: Resistivity Vs Doping concentration for TE, and Experimental Values, L=1220A° 



Dopant concentration (cm' 3 ) 

Figure 14b: Resistivity Vs Doping concentration for TFES+TE, TE and Experimental Values 




Numerical Calculations And Results 


56 


S.No 

Dopant 

Concentration 

[cm -3 ] 

For Grain Size= 12204° 

Experimental 
Resistivity [flcm][ 
7,20] 

TFES+TE 

Resistivity 

Only [Q.cm] 

TFES 

Resistivity 

[Q.cm] 

TE 

Resistivity 

[Q.cm] 



EHEH 



7.2400X10^ 


3.9811X1016 

1. 2589X1 0 US 



2.1 400X1 0 U5 


1.1749X1017 

5.2481X1 0 W 

4.4668X1 0 W 

3.4700X1 0 °° 

IrTTrnriM 




5.4954X1 0 U1 

4.9000X10^ 

HRSIiliFXIiyBl 


3.1261X1017 

1.81 97X1 0 UI 

1. 3804X1 0 UI 

7.5858X1 0 U1 

MU-liMliH 


9.2683X1017 

2.4547X1 0' u1 

9.5499X1 0^ 




8.9125X1018 

2.5704X1 0' 1 * 

5.8884X1 O^ 

2.3442X1 0^ 

1.1220X10 


5.0000X1019 

1.2303X10'^ 

2.6303X1 O' 00 

1. 0233X1 0"^" 

4.6774X1 0' UJ 


Table 12: Concentration and Resistivities due to different effects. 
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4.2.8 Plot and Corresponding Data of “Resistivity vs. Doping 
Concentration” for Different Grain Sizes (200A°-1 220A°) 



Figure 15a: Resistivity Vs doping concentration for different grain size 
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Table 13: Concentration and Resistivities for Different Grain Size 
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4.2.10 Plot and Corresponding Data of “Current Density Vs Voltage (J-V) 
Plot” 


Current density Vs voltage characteristic is drawn using 
their combined relation given equation (43a), 

( V / 2 


j = 


2 

q N 


2rnnfjkT 


-qV b /kT 

e V 

o 


both TFES and TE, 


i.e. 


V 0 is the applied voltage across the composite grain boundary region. If the 
applied voltage is V a and the number of grains are g. 



J = q N, 


ImrifokT 


1/2 


-qv b lkT 


V a 

g 



Figure 16: Current Density Vs. Voltage at Room Temperature 
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S.No. 

Voltage ( V) 

Current Density (Amp.cm^) x10 / 

1 . 

1 

.0407 

2. 

2 

0.0815 

3. 

3 

0.1222 

4. 

4 

0.1629 

5. 

5 

0.2036 

6. 

6 

0.2444 

7. 

7 

0.2851 

8. 

8 

0.3258 

9. 

CD 

0.3665 

10. 

10 

0.4073 

11. 

20 

0.8145 

12. 

30 

1.2218 

13. 

40 

1.6291 

14. 

50 

2.0363 


60 

2.4436 

16. 

70 

2.8509 


Table 15: Voltage and Current Density 



























Chapter 5 
Analysis And Discussion 


5.1 Discussion About Grain Boundary 

The theoretical model was matched to experimental data (Table 1-4) using the 
width of the grain -boundary barrier, W gb , as fitting parameter W gb \s determined 
by using Newton-Raphson iterative method. The value of W gb comes out to be 
1 1 .86 A 0 at trapping density 3.34x10 12 cm' 2 (Table 1) and at very high doping 
concentration it comes out to be 1 1 .89 A 0 (Table 3). It should be noted here that 
the value for the W gb is physically reasonable. A grain-boundary height of 0.66 
eV gives a grain boundary band gap of 1.32eV, which satisfies the basic 
assumption of the model. A grain boundary thickness of 11.86A 0 or 11.89 A 0 is 
roughly the width of two -three mono layers, which seems to be reasonable 
when compared to grain boundary width of 7 A 0 in case of phosphorous doped 
polysilicon[4]. Grain boundary is the transitional region separating the highly 
ordered crystallites. The grain boundary width is different for P, & As [4] as 
compared to B, hence it can be concluded that grain boundary width is the 
property of dopant species. It is also seen that, both at moderate and high 
doping concentration there is not a significant change in the W gbi implying that it 
is not a function of dopant concentration as well. 

5.2 Discussion of Various plots 

1. In section 4.2.1, Figure (7), plot of barrier height vs. doping concentration 
is shown; the barrier height increases with doping concentration and after 
' reaching a maximum value starts decreasing. It is due to the fact that 
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barrier height varies differently with doping concentration at low and high 
doping concentrations respectively. When the doping concentration is 
lower than the trapping density ( LN<Q t ), the crystallite is completely 
depleted of carrier and the traps are partially filled. According to Seto [2] 
Barrier height (V b ) increases linearly with doping concentration. Linear 
variation is given by, 



where N is the doping concentration. 

Seto [2] assumed that grains are completely depleted at low doping 
concentration, which is possible only at very low doping concentration not 
at moderate doping. Therefore the relation (2W=L) distorts the linear 
variation at low doping concentration as shown in plot of Figure (7) 

When all the traps are filled ( LN>Q t ), only part of the crystallite is depleted 
of carriers. V b varies inversely with doping concentration, thereby 
decreasing the V b at high doping concentrations. 

b 8 sN 

Our plot shows the same potential barrier values as shown in [21,22] 
experimentally. 

Figure (8a) shows the plot between Mobility Vs doping concentration. 
Since the Mobility depends on the negative exponential of V b (equation 
43c), therefore, when the V b increases mobility decreases. Our results are 
similar to that in [2,21]. Value of mobility at very low doping concentration 
is shown in Figure (8a), we see that curve is almost flat, that means at 
very low doping concentration there is not appreciable change in mobility. 
Figure (8b, 8c) shows that contribution due to thermionic emission deviate 
from the effective mobility value at high dopant concentration. These plots 
also show that effective mobility value reaches to single crystal range at 
high dopant concentration. The reason is clear from equation (39) as 
thermionic mobility is very high so its reciprocal is negligible compare to 
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single crystal. Therefore, mobility at high dopant concentration is 
approximately in single crystal range. This result is same as given in [21]. 
In sections 4.2.2 Figure (9a, 9b) and in section 4.2.3 Figure (10a, 10b) 
plots of mobility Vs doping concentration are shown and a comparison is 
shown between (TFES+TE), TFES, TE and the experimental values for 
two different grain sizes 200A° and 1220A° respectively. We analyze that 
at low doping concentration thermionic emission (TE) matches with 
experimental values. In section 3.3.2 reason is given why at low doping 
concentration TE dominates. At high doping concentration (10 17 -10 18 ) 
contribution of (TFES+TE) matches with the experiments (Figure 9-10). 
But as the concentration increases (1 0 18 -1 0 20 ) contribution due to TFES 
increases significantly The reason is clear from the fact that all our 
assumptions are satisfied at high dopant concentration. Similar results are 
obtained for both low and high grain sizes. Better results are obtained in 
case of higher grain size. Our results predict the same conclusion as 
given in [2,7,20,21,24]. 

2. In section 4.2.4, Figure (11), plot of mobility as a function of (VkT) is 
shown, based on empirical equation (43f), which shows a straight-line 
behavior with negative slope equal to the barrier height (qV b ). Plot is 
drawn at barrier height 0.078V, grain size 1220A° and doping 
concentration 8.92x1 0 17 crrf 3 . Figure (11) is provided for making the 
comparison with the values obtained from the empirical relation. 

3. In section 4.2.5, Figure (12a, 12b), shows the plots between doping 
concentration Vs active carrier concentration. We know that at low doping 
concentration majority of the carriers are trapped in the grain boundary; 
hence the active carrier concentration is very small in comparison to the 
doping concentration. When all the traps are filled majority of the carriers 
contribute to the conduction. Therefore active carrier concentration 
increases at high dopant concentration as shown in the Figure (12a). At 
very high doping concentration active carrier concentration is almost same 
as doping concentration. We see in the Figure (12a, 12b), that there is a 
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increase in active carrier concentration between doping range ~10 17 - 
2x1 0 18 . This is result of the fact that as traps gets nearly filled more & 
more dopant atoms contribute to the free carriers. Further increase in 
dopant concentration reduces the depletion layer and eventually the free 
carrier contribution approaches that of crystalline silicon this is also in 
agreement with results in [2] and [21]. 

4. In section 4.2.6, Figures (13a, 13b), and in section 4.2.7, Figures (14a, 
14b), show the plots between resistivity Vs doping concentration for low 
as well as high grain sizes. The contribution due to TE, TFES, (TFES+TE) 
is compared to the experimental values. At low doping concentration plot 
of TE alone matches with the experimental values. At the high doping 
concentration (10 17 -10 18 ) contribution of (TFES+TE) matches with the 
experimental values, while TE considered alone deviates. But as the 
concentration is further increased from (1 0 18 -1 0 20 ) contribution of TFES 
increases significantly. Therefore in this range contribution of TFES alone 
matches with Experimental Values (Figure 13-14). We got similar results 
as shown in [2,7,20]. Thus we see that resistivity cannot be explained 
using a single mechanism of transport. It is governed by different 
phenomenon at low and high doping concentration. These results give the 
conclusion that TE is dominant phenomenon at low doping concentration 
and TFES is a dominating phenomenon at very high doping 
concentration. 

We see in Figures (13-14) that at low doping concentration resistivity is 
very high then an abrupt change occurs in resistivity between doping 
range 10 17 -10 18 (L=1220A°) and 10 18 -10 19 (L=200A°) after that resistivity 
decreases approximately linearly. The reason of sharp decrease in 
resistivity is clear from the discussion above that under this range number 
of free carrier increases abruptly, consequently resistivity decreases .The 
reason of high resistivity at low doping concentration is that, most of the 
carriers are trapped leaving only a few free to contribute to the 
conduction. Carrier trapping results in the reduction in the number of free 
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carriers within the grains. Hence resistivity is very high at low doping 
concentration At high doping concentration (10 18 -10 19 ) the resistivity 
reaches to single crystal range as illustrated by Table 5. 

In section 4.2.8, Figure (15a), plots of resistivity Vs doping concentration 
for grain sizes 200A°and 1220A° are shown. We get the lower values of 
resistivity for higher grain size. Similarly in section 4.2.9 Figure (15b) plot 
of resistivity Vs different grain size is shown, resistivity decreases with 
increasing gain size. Reason of lower value for higher grain size is 
obvious from equation (42b, 43d), resistivity is inversely proportional to the 
grain size. 

5. In section 4.2.10 Figure (16) plot of current density Vs voltage is shown 
using equation (43a), Numbers of grains are taken 100 from Baccarani 
[22]. We got a similar hyperbolic sine J-V curve as shown in [5,20,22,]. 
Equation (43a) can be used to calculate the number of grain. 



Chapter 6 
Conclusions 


This work was carried out to study the electronic properties of boron doped 
polysilicon films. This material has drawn considerable attention for its potential 
application in polysilicon devices namely n-channel thin films transistors (TFT s). 
The objective of the studies was to have better understanding of electronic 
properties of the material namely mobility, conductivity (resistivity), transport 
mechanism in general and how the free carrier concentration varies as a function 
of doping levels in the material. 

The study concludes that for fine grain Boron doped polysilicon (grain size 200- 
1220A°), at high doping concentration (>10 18 cm 3 ), most of the dopants 
contribute to free carriers. However, at lower doping concentration only a fraction 
of the dopants contribute to average carriers. This result, though previously 
reported by other researchers up to 10 16 cm 3 , did not cover the doping range as 
low as 10 14 cm' 3 . The study shows that at doping concentration of 10 14 cm 3 , the 
free carrier concentration is ~10 9 cm 3 . Similar results were obtained for 200A 
and 1220A° grain sizes. 

The study also concludes that et low dopant concentration (<10 cm ) 
Thermionic emission (TE) is the dominant mechanism behind carrier transport 
across the grain boundaries. However, both thermionic emission (TE) and 
thermionic field emission with scattering (TFES) contribute to carrier transport 
across the grain boundaries at high doping concentration (10 -10 cm ). 
Further, at very high doping concentration TFES is the dominant mechanism of 
transport across the grain boundaries. 

The study also shows that resistivity decreases with increasing grain size, which 
is in general agreement with the works of others. 
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Scope for the future work 

In the present work assumption has been made that the applied voltages is less 
than kT/q, so that in addition to thermionic emission carriers tunnel through the 
grain boundaries. However, at high voltage TFE may also contribute to the 
conduction mechanism and must be considered. 



Appendix A 

Potential Barrier Equation 


Writing the Poisson’s equation for depleted grain 

J 2 F(x) _ qN 
dx 2 s 

On integration, 

?m=*jL x+Ci 

dx £ 


W dV(x) 

After applying the boundary condition at x-W + — = 0 as the field is 

zero in the neutral region. 


dV(x) 

_qN 
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( w Y] 
w 

dx 
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Again integrating we get, 


V(x) = 


qN_ 

l£ 
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Now at x = —^-,V(x) = V b 
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( w \ 
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l 2 J. 
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.At x=0 the first term on the right hand side of the above equation becomes zero 
and potential is also zero. 
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At low doping concentration grains are completely depleted, so 2W = L. 


n = 


qNL 2 

Ss 


Appendix B 
W.K.B Method 


This method applies only situations in which the potential is slowly varying 
function of position. 


Let ¥(*) be the wave function satisfying Schroedinger's equation. 


2m 


8x 


r 




(i) 


again, f the solution of equation (1) be in the from ¥ = Ce^ x),h , then we get from 
equation (1). 


ih<f>" (x)~ <j>' 2 (x)+ 2m [E - qV (x)] = 0 


( 2 ) 


To get the approximate solution of (2), we apply W.K.B method and hence 
expand ${%) in powers of ft i.e. 



as fl is very small , the first two term in the above equation give sufficiently 
good approximation of • So putting the value of in equation (2) and 


collecting coefficient of various power of ft . We get , 


t t 


2 m[(E~qV(x))-</>Q] + h i<t>o ~~ 2( f*oh 


+h J 


it *2 1 1 

i<f 1 -<f\ 


= 0 


( 4 ) 
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thus from equation (4) we get , 


0 o =±J[2m{E-qV(x) 


:.<p 0 =±lj\2m{ 


<t>\ = + c l 


\Sx ) 1 b , m ' dx 

2 2[2 m{E-qV(x))f 2 4 a [2m(E-qV(x)f 2 


, M SV 

We see that yiy^) will be small whenever ~ is small and [E-qV(x)\ is not 


too close to zero. 


^W«^oW+y ln K 

r . b ^ ? ^ 

v F(x) = C[2m(E -qV (x))] -1/4 exp ± ^ J ^[2m(E-qV(x))\dx ^ 


This ± show the direction of motion of the wave function in positive & negative 
direction. 

We can write it in more compact from using the momentum expression. 


p = y/2 m\ 


C r Up 


v i , (x)= -p=-exp ±i]ydx 
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Transmision Pr obability(D)= 


Transmittedlntensity 

Incidentlntensity 


After solving transmitted intensity & incident intensity we get, 


|c| 

Transmitte dlntensity = L ~ L - 

P 


Incidentln tensity 


|cf 


Pi.. I 1 


e Xp j,_ dt |._ exp |,_ A j 


■\ 


i 2 y (8) 
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Under W.K.B approximation, j h 


Thus; 


( 9 ) 




Appendix-C 
Concentration Density 


Say 


*(£,) = 

h m, o o 
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1 + exp 


' 2 2 
P x + Pp 
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Solving the first integral. 
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Now for the tunneling to be occurring K.E of the holes should be lesser than 
Fermi energy. 

2 

— "V + s < E F So with this condition and 
2 m h 
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\ 2m h 
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After solving & putting the proper limits, 
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Appendix-D 
Taylor Series Expansion 


Taylor series gives you the expansion of the given function around to a close 
point. 


f(x + a)= f{x)+af\x)+—f"(x)+... 


We have tunneling probability. 


D{E X ) = exp 


or 


- 4 n b r 


f[ 2 ™* h {qV (x)-E x )ldx\ 


In D(E X ) = 


- 4/r * L * 


^ r « a 

)\lm* h {qV {x)~ E dx\ 


(1) 


So using Taylor series we can write, expand the R.H.S of equation (1) around 
the energy s x =qV b -E x so we have. 

- An 


f(E x ) = | -^-j[2ml(qV(x)-E x )]2dx 


Or f(qV b -s x ) = j[2 m' h (qV(x) - (c qV b - s x ))]* & | (2) 

But according to Taylor series expansion equation (2) is given by , 

AqV, -e„) = fm-fm'e, + /(?r»)V - (3) 
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Value of /(#F 4 )can be calculated by putting q-V b qV b +s x on both side of 
equation (2), 

/(^a) = j— fym^qVW-qV^dx 

To get f(qV b ) differentiating equation (2) w.r.t. qV b & Again putting 
qV b -> qV b +s x on both side of differentiation. 

f(qV„ -£,)'= )[2 ml(qVU) - ( qV t 

f(qV t ) =\^f)[2ml{qV(x)-qV b )fdx 

Similarly the second derivative can be calculated 

C b -3 

fiqvS = ^ l[2ml(qV(x)-qV t )]Tdx 
Now the equation (3) can be written as follows, 

f{q Vb ~ £ x ) = — IT 1 ! + C \ £ x + f\ £ x ] 

Where bi, Ci & fi are given by, 

6, m , i { 3 V(x)-qV b ^ix 

c l =\^-)[2m' i (qV(x)-qy l ,)fdx 

/, =k-){2m’ h (qr(x)-qr t )fidx 

So finally equation (1) is given by, 

In D(E X ) = -/■&?+ Ci£ x + fis x 2 + 7 

or D(E x )=exp-[ bi+ c^Sx + Us 2 + 7 
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Value of constants 



We have to solve b\ & Ci for carriers in the energy range qV b <E x <qi/>, so taking 
the value of qV(x) from equation (13) we get . 


qV(x) = 0 
qV(x) = q</> 


for H>r„/2 
for H<^/2 


-'r 


>- 




Which after solving and applying the limits gives 

b x =^-W sb [2ml} l2 {q l f-qV b )'l 2 
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Similarly the value of the c ? can be calculated, 




so using this integral & putting the same limit we get the value of C\. 




Appendix-F 
Current Density 
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Appendix-G 
Potential Across Grain Boundary 


Writing the Poisson’s Equation for the grain boundary region, we get, 


ForH< 5t 

dx 2 e 11 2 


( 1 ) 


w 

Integrating and applying the boundary condition at x = , V(x) = <j) 
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dx s 


< ^ 

x — 
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Further integrating and applying the boundary condition at x = = (j> 
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2s 
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As V d is half of the potential drop across the depletion region. It can be 
calculated by taking the difference of above two equations. 
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*K N G 


4s 


or 


K= £t (ase ‘= w « N ° : 


(4) 


Now from the [23] dealing our grain boundary likes metal —semiconductor 
(Schottky barrier) problem. We can write the expression of total potential across 
semiconductor & at the interface of grain boundary. 


V = K + V, 


v*\ + WgbqD * s v + 


qNc 

\2s j 
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(5a) 

(5b) 


where D ss is the hole diffusivity in the silicon grain 

Equation (5a) & (5b) represent the same potential ,so the potential across the 
grain boundary can be found by taking the difference of last term in the 
expressions (5a)& (5b) ,i.e. difference in potential at interface. 


y *~ 


f M 


qNo 
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( 6 ) 


Now as the applied voltage is very less as compare to barrier height (F 0 « V b ) 
V L = V b And V R = 2V d - V b by [4], with these relation equation (6) changes to as 
follows. , 


.W.[€o 

V 2s 




v gb =w. 


(2^) 
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So finally equation (4) and (7) give relation between potential across grain 
boundary and total potential drop across the depletion region. 


V^=W. 
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Appendix-H 
Newton-Raphson Method 


Let x 0 be an approximate value of a root of the equation f(x)= 0, and let x 0 +h be 
the exact value of the corresponding root, where h is very small quantity. Then, 

f(xo+h)= 0 (1) 

Since xo+h is the root of the equation f(x)= 0. 

Expanding (1) by Taylor’s theorem, we get 
f(x 0 +h) = f(x 0 ) + hf'(xo) + h 2 f '(x 0 )/2! + = 0 

Since h is very small, neglecting second and higher order terms and taking the 

first approximation, we have 

f(x 0 ) +hf (x 0 ) = 0 => h =-f(x 0 )/f (x 0 ) , provided f (x 0 ) * 0. 

.-.x^ x 0 +h = x 0 - f(x 0 )/f (x 0 ) (2) 

Relation (2) gives the improved value of the root over the previous one. Now 
substituting Xi for xo and x 2 for xi in (2), we get 

x 2 = xi-f(xi)/f(xi) (3) 


In general, we can get an approximation or the iteration formula 


X n+ 1 = X n - f(X n )/f(X n ) 


( 4 ) 


From this formula, we can calculate successive better values of the root. 
Formula (4) is known as Newton-Raphson Method. 
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